Common path interference in Zener tunneling is a universal phenomenon 
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We show that the probability of electric field induced interband tunneling in solid state systems 
is generically a non-monotonic (oscillatory) function of the applied field. This unexpected behavior 
can be understood as arising due to a common path interference between two distinct tunneling 
solutions. The phenomenon is insensitive to magnetic field, and arises whenever the low energy 
dispersion relation contains higher order terms in addition to the usual p 2 term. Such higher order 
terms are generically present, albeit with small co-efficient, so that the oscillatory Zener tunneling is 
a universal phenomenon. However, the first 'Zener oscillation' occurs at a transmission probability 
which is exponentially small when the co-efficient of the higher order terms is small. This explains 
why this oscillatory aspect of Zener tunneling has been hitherto overlooked, despite its universality. 
The common path interference is also destroyed by the presence of odd powers of p in the low energy 
dispersion relation. Since odd powers of p are strictly absent only when the tunneling barrier lies 
along an axis of mirror symmetry, it follows that the robustness of the oscillatory behavior depends 
on the orientation of the tunneling barrier. Bilayer graphene is identified as a particularly good 
material for observation of common path interference, due to its unusual nearly isotropic dispersion 
relation, where the p 4 term makes the leading contribution. 

PACS numbers: 



I. INTRODUCTION 

Tunneling and interference are two signature quantum 
effects that are commonly thought to occur in very differ- 
ent experimental setups. Interference usually involves the 
superposition of distinct saddle point solutions, whereas 
tunneling through a barrier proceeds through the evanes- 
cent solution with maximum decay length. Recently, it 
was shown that particles in bilayer graphene exhibit in- 
terference underneath a tunnel barrier, leading to a Zener 
tunneling probability that is a non-monotonic function of 
barrier height [1] . The interference was further identified 
to be of a common path type, and hence unaffected by 
magnetic fields. The tunneling intensity is an oscillatory 
function of the bandgap and the potential: 

T = 4|a| 2 eT^ CO s 2 (jr + <t>) W 

where S r and Si are monotonic functions of the bandgap 
and field strength, and a and <f> are constants. The 
bandgap and field strength can be tuned so that the tun- 
neling current becomes zero. 

In this work we show that such common path inter- 
ference underneath the tunnel barrier is not specific to 
bilayer graphene, but is in fact generic, occurring in any 
system where the dispersion relation contains a p A term 
in addition to the usual p 2 term. Since higher deriva- 
tive terms generically arise in effective low energy Hamil- 
tonians 5 , the common path interference is universal. 
However, if the coefficient of the high derivative term is 
small, then the first 'Zener oscillation' does not manifest 
itself until an exponentially small tunneling probability is 
reached, explaining why this common path interference 
effect has hitherto been overlooked. The unique feature 



of bilayer graphene is that the higher derivative terms do 
not come with small coefficients, and thus the Zener oscil- 
lation phenomenon should be (relatively) easy to see. We 
also show that there is nothing special about a p 4 term, 
and that common path interference generally arises when 
the dispersion relation contains higher powers of p 2 (but 
the odd powers of p are absent). Common path inter- 
ference may be detected either in transport experiments, 
or in optical absorption experiments designed to detect 
photon assisted tunneling. Bilayer (and ABC) multilayer 
graphene are identified as particularly promising candi- 
dates for detection of this phenomenon, but other narrow 
gap semiconductors (such as several recently discovered 
topological insulators) may also offer a potential experi- 
mental playground. 



This paper is structured as follows. We begin in sec- 
tion II by explaining the origin of the Zener oscillations 
using the WKB approximation. We generalize the anal- 
ysis in pP to a dispersion relation with a quadratic piece, 
and demonstrate that the common path interference sur- 
vives. In section |III[ we show that in fact common path 
interference is a generic feature for dispersion relations 
that contain higher powers of the momentum. Finally, 
in section III, we verify our results (for the fourth order 
polynomial dispersion relation) by re-deriving them in an 
independent manner, by mapping the problem to a differ- 
ential equation which can be numerically solved. Finally, 
in section IV, we present our conclusions and comment 
on the possibilities for observation in experiment. 
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II. WKB APPROXIMATION 

In this section we calculate the probability of Zener 
tunneling as a function of applied field using the Wentzel- 
Kramcrs-Brillouin (WKB) approximation. A convenient 
toy model to investigate the physics we have in mind is 
the two band effective Hamiltonian for bilayer graphene 
013! 
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Here, A = E g /2, where E g is the band gap, which can 
be tuned using external gates. The work [T] analyzed 
the above Hamiltonian at = 0, however, a rigorous 
treatment of the effects of higher bands generates a small 
but nonzero j3 < PJ. Note that (3 < leads to bands 
that have a 'mexican hat' shape, whereas f3 > leads to 
bands that have only one minimum. A great advantage of 
the Hamiltonian above is that the low energy dispersion 
has the form 



\E\ = A 



(2) 



which is the general polynomial dispersion of interest to 
us in this paper. A quadratic dispersion (a„>2 = 0) is 
conventional, and arises for e.g. free fermions governed by 
the Schrodinger equation. The tunneling behavior in this 
limit is well understood, and does not display common 
path interference. However, terms at higher order in the 
momentum are generically present in the dispersion rela- 
tion about the bottom of any bandstructure (other than 
free fermions) as can be explicitly seen, e.g. by perform- 
ing k.p perturbation theory for order higher than two 
[5] . This paper is about what happens when these higher 
momentum terms are nonzero. It is sufficient to include 
just the first non-trivial term, ai ^ 0. We will demon- 
strate that this leads to common path interference, and 
a tunneling probability that is a non-monotonic function 
of barrier height. Including further terms 03,04... does 
not change the physics, as we show in Section III. We 
note that an added advantage of the Hamiltonian (|TTJ) is 
that ai oc j3. Thus, by taking /3 — > 0, we can recover the 
model analyzed in lj, which serves as a useful double 
check on our results. 

We will consider the Hamiltonian (|H| with a one- 
dimensional linear potential V(x) = Fx, where F is the 
force on an electron due to an (in-plane) electric field. 
The treatment here can be easily generalized to other 
monotonic potentials. In the presence of this potential, 
Zener tunneling of particles between conduction and va- 
lence bands will occur [3]. We calculate the probability 
of Zener tunneling using a WKB approximation, in the 
manner of [1]. 

For a slowly varying potential along the s-axis, the 
WKB approximation states that the change in the wave- 
function after propagating from xq to x is proportional 



to e 1 ^o Px ^ dx . To find the amplitude of the tunneling 
through the barrier, we have to solve for p x and integrate 
it over the forbidden region of the energy gap, in which 
it has a non-zero imaginary component. 

For simplicity let us consider just electrons at normal 
incidence, p y — 0. The inclusion of a non-zero p y com- 
plicates the analysis without changing the physics. We 
will focus on electrons at normal incidence in this pa- 
per. Then the momentum along the x-direction, p x = p. 
Finding eigenvalues of H + V(x) yields the equation : 



P 



2A/3 



1 

iiu- 



;P 



1 

4 m- 



■0 2 







(3) 



where E' = E — V(x). Let us define for brevity the 
quantities 



1 02 A P r A 2 - E' 2 

7=-, S=— . (4) 
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Therefore, p 4 + 2"fp 2 + 5 = 0. The solutions to this 
quadratic equation in p 2 are: 



P 



= -7 ± \/"f 2 - 5 



(5) 



This equation implies that there are three distinct re- 
gions with sharply different physical properties (Fig. 1), 
and we will discuss each region in turn. The essential 
ideas are illustrated in Fig.l (for /3 < 0) and Fig. 2 (for 
(3 > 0) respectively. 



A. Region I: 5 < 

In this region there are 2 real solutions: 



p x = ±y -7 + \ll 2 - & 
and two imaginary ones: 

Px = ±iy 7 + Vl 2 - S 



(G) 



(7) 



The former two solutions correspond to propagating 
waves moving left and right respectively, while the latter 
correspond to evanescent waves. This region corresponds 
to the region outside the tunneling barrier. 



B. 



Region II: < 5 < 7 



For j3 < 0, this region is outside the tunneling barrier. 
However, there are four real solutions given by 



Px = ± 



-7 



±W 2 



(8) 



To understand the existence of four solutions, we have 
only to remember that the bandstructure for /? < is 
a mexican hat. A line cut at low energy will intersect 
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the mexican hat at four points, and these four points 
correspond to the four real solutions identified above. 

Meanwhile for f3 > 0, there are 4 imaginary solutions 
given by 



r 



(9) 



Thus, for f3 > 0, this region is inside the tunnel barrier. 
However, in this region the tunneling solutions are non- 
degenerate (have different decay lengths), and tunnel- 
ing will be dominated by the solution with largest decay 
length. 



C. Region III: S > j 2 

In this region, the solutions are all complex, and take 
the form 



Px 



= ±(~7±«V / -7 2 + '5) 1/2 
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(10) 



The lack of real solutions means that this region is inside 
the tunneling barrier. Only two of the solutions above 
correspond to tunneling from left to right. However, both 
are degenerate (have the same decay length). As a result, 
both will contribute equally to tunneling, giving rise to 
common path interference. 



D. WKB solution 

The real and imaginary parts of p are depicted pictori- 
ally in Figs. [l]and[2]for a negative and positive value of f3 
respectively. For tunneling from left to right in these fig- 
ures, we choose solutions so that p x is always continuous. 
Also, in regions II and III, any imaginary component of 
p x should have positive sign, so that that the solution de- 
cays to the right. Under these constraints, the possible 
solutions for j3 < (the relevant case for BLG) reduce to 

Region I: p x = ±^-7 + \/"f 2 -5 
Region II: p x = ±y-7 + \Jl 2 - 5 



Region III: p x = ±yj + ^ 
Note that in region III there are two degenerate tunneling 
solutions. We will show that interference of these solu- 
tions will lead to an expression of the form ([!]). Mean- 
while, for (3 > we obtain the solutions 

Region I: p x = ±\J -j + \/j 2 -S 



Region II: p x = i\J -7 ± ^j 2 - 5 
Region III: p. 
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Again there are degenerate tunneling solutions in region 
III, which will give rise to common path interference. 
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FIG. 1: Real (above) and imaginary (below) parts of p x at 
2m/? — — l.Red circles are \J — 7 + \Ay 2 — S, green squares 
are 



- \J — 7 + sj^j 2 — <T, blue diamonds are 
and purple stars are — y— 7 — \/7 2 — 



We now explicitly evaluate the tunneling probability. 
The end-points xl and xr of the barrier region are the 
points at which 7 2 = 8. This yields: 
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(11) 



Note that here any form of V(x) can be used to obtain 
the endpoints, as long as it is smooth and monotonic. 

For /3 < 0, the barrier is only in region III. The change 
in the wavefunction after it has moved through region III 
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FIG. 2: Real (above) and imaginary (below) parts of p x at 
2mj3 = 1. Red circles are y— 7 + x/'y 2 — 5, green squares are 

-y — 7 + \/7 2 — ^, blue diamonds are y— 7 — \/7 2 — 5, and 



purple stars are —v—l— \J l 2 — & 



is proportional to e /3b (e l/3a + e~ l/3a ), where J 3a and 7 3 b 
are the integrals over the barrier region of the real and 
imaginary part of the momentum respectively. We now 
explicitly evaluate these integrals. We have 
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A 2 - E' 2 A/3 



dx (12) 



Changing variables to E' = E + Fx, we obtain 
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(13) 



where E' L = E + Fx L = ^= and E' R = E + Fx R = 
2m^/c ' After scaling out all dimensionful quantities, this 
can be recast as 

2A 3/2 m l/2 
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and On 



= 

cos~ 1 (|2m/3|/-y/l + (2rn/3) 2 ) are the two fixed pa- 
rameters that characterize a particular tunneling 
process. The function a 3a is plotted in Fig. 3. 
Similarly, is given by 
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Again, a 3 fc is plotted in Fig. 3b. Note that a 3a and a 3 b 
are both strictly positive. 

Meanwhile, for j3 > 0, for the tunneling through region 
III will have exactly the same form as for j3 < 0, except 
that a 3a and will interchange with each other, i.e. 
«3a(/? > 0) = a 3fc (/3 < 0) and a 36 (/3 > 0) = a 3a (^ < 
0). Also there will be a contribution to the decay of 
the wavefunction amplitude as it moves through region 
II which is of the form 
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where a 2 is plotted in Fig. 3. After tunneling, the wave- 
function will be multiplied by a factor e ~( l3a+2l2 \e ll3b + 
e~ ll3b ). Note that once again we will have common path 
interference coming from region III, which will give rise 
to an oscillatory tunneling probability ([I]). 

Thus, we have shown that common path interference 
occurs and results in tunneling oscillations whenever the 
dispersion relation has a quartic component. The exis- 
tence of common path interference does not depend on 
the sign of the quadratic component. However, as the 
ratio of quadratic to quartic co-efficients is made big- 
ger, the tunneling oscillations are pushed to larger and 
larger bandgap, and the first antinode occurs at exponen- 
tially smaller tunneling amplitude (Fig. 4). This follows 
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FIG. 3: ai, Q2 and 03 vs |2m/3|. 



because the exponential suppression of tunneling is pro- 
portional to 0:3, and a.3 goes to zero as 2m /3 — > 00. 

In most materials, the co-efficient of the quartic term 
is much smaller than the corresponding coefficient of the 
quadratic term. As a result, the Zencr oscillations, while 
theoretically present, will be extremely weak, and may 
not be visible experimentally. The special feature of bi- 
layer graphene is that 2mf3 <C 1, and as a result the Zener 
oscillations should be relatively easy to see. 



III. GENERALITY OF COMMON PATH 
INTERFERENCE 

In this section we demonstrate that common path in- 
terference is not simply a property of quartic dispersions, 
but in fact is a general feature for any dispersion of the 
form 
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(17) 



with real coefficients a„ 
proceeds as follows. 



The argument is simple and 



The equation (17) is a N th order polynomial in p 2 
with real co-efficients. By the complex conjugate root 
theorem, complex solutions must come in conjugate 
pairs. Additionally, in the barrier region, there cannot 
be a real positive solution for p 2 . This implies that 
for even N, the solutions to (17) must have the form 
P 2 = {rje ±l2 ^ j }, where j = 1...N/2 and r and are 
real numbers. The corresponding solutions for p are 



ing to tunneling from left to right are p — { ^Jrj (± cos c^j t 
% sin <f>j)}, which all come in degenerate pairs. The system 
simply picks out the degenerate pair with longest pene- 
tration length (minimum ^/rj sin <f>j ) . Since there are two 
solutions with this decay length, the two solutions can in- 
terfere, to give rise to the physics discussed above. Note 
that this is true for any (even order) polynomial disper- 
sion relation, not just the quartic dispersion relation we 
discussed thus far. 

For odd N the story is similar, except that there is also 
one solitary solution p 2 = —A 2 , where A is a real num- 
ber. This leads to solutions p = ±iA. Common path 
interference can be absent for odd N if and only if A 
is the longest penetration length in the problem i.e. IFF 
A < min ^rj sin . If this very special condition is satis- 
fied then there will be no common path interference, and 
this single non-degenerate mode will dominate tunneling. 
Otherwise the story will be similar to that for even N, 
and there will be common path interference. 

Thus, we have shown that for a dispersion relation that 
is an N th order polynomial in p 2 (with real co-efficients) , 
there must necessarily be common path interference if 
N is even. If TV is odd the system may escape common 
path interference, but only if the single non-degenerate 
tunneling mode has longer decay length than the (N — 1) 
degenerate tunneling modes. 

We have focused on electrons at normal incidence, 
since tunneling at oblique angles is exponentially su- 
pressed by the longer tunneling path. We have also as- 
sumed that the dispersion relation does not contain odd 
powers of p (i.e. the dispersion relation is symmetric un- 
der p — » — p). This is rigorously true for a tunneling 
barrier parallel to an axis of reflection symmetry of the 
underlying crystal. 

For a general tunneling barrier, however, odd powers of 
p will arise in the dispersion relation due to crystal fields, 
albeit with small co-efficient. These odd powers will lift 
the degeneracy of the tunneling solutions, and will ulti- 
mately destroy common path interference. However, if 
the co-efficients of the odd p terms are small, then the 
degeneracy will be barely lifted, and the degeneracy lift- 
ing will only manifest itself on lengthscales much larger 
than the tunneling decay length. In this event, common 
path interference should still be visible if the barrier is no 
too wide - i.e. we should still be able to see the first few 
oscillations, but high order nodes will be washed out by 
the lifting of the degeneracy of the two interfering modes. 



IV. SOLUTION BY MAPPING TO 
DIFFERENTIAL EQUATION 

An independent way to verify the conclusions from the 
WKB treatment in real space is to consider the prob- 
lem in momentum space. Following Ref. [T], we substi- 
tute x = ihdp x . This allows us to write the Schrodinger 
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FIG. 4: Transmission as a function of the bandgap for differ- 
ent values of 2m/3, E g measured in units of A . The points 
are numerical data and the lines are fits to the points with 
functions of the form [I] S r and Si are obtained from the 
WKB approximation while a and <f> are chosen for the best 
fit. 



equation Hty = Eip as 
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where the <jj are the Pauli matrices, and we have taken 
E = by choosing the origin of coordinates appropri- 
ately. This equation can be interpreted as a two state 
system swept through an avoided crossing, with p x play- 
ing the role of time. The probability of Zener tunneling 
is equal to the probability that a system initialized in a 
particular eigenstate ends up in the orthogonal eigenstate 
after the sweep through the avoided crossing. Again, we 
focus on electrons at normal incidence, p y = 0. 

For |2to/3| < 1, as p x = ±oo, the eigenstates of the 
Hamiltonian are eigenstates of cr\, i.e. <\>\ — [1,1] and 
<t>2 = [1,-1]. The Zener tunneling probability is equal to 
the probability of starting in the state and ending in 
the state </> 2 after the sweep through the avoided cross- 
ing. In a simple approximation, we can argue that the 
transition has a chance of happening when the coefficient 
of the <73 term is larger than that of the a\ term: 



< A + /3p 2 x 




PA 



(19) 



C202- | Ci | 2 + | C2 1 2 = 1. This gives us the two simulta- 
neous differential equations: 



dp x 



ihF 



dC 2 
dp x 
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(20) 
(21) 



with the initial conditions at p x = — pa, |Cb| = 1 and 
|d|=0. 
These are solved by 



C-2 



i cos I ui I p x + -^r-p x I + A 
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(22) 
(23) 



where, uj = From initial conditions, A 

UPA (l + 3^Pa) • 

The total phase acquired at p x — p\ is 



A, 



3/2 



a{m/3) (24) 



Let the wavefunction during the transtion be ^ — C\4>\ 



which indicates that the probability of tunneling ( | C 2 1 2 ) 
will oscillate with period dependent on m/3. 

For |2m/3| > 1 however, the coefficient of the 03 term 
is always larger than that of the <j\ term. The simple 
approximation above is no longer applicable then since 
at no time does the transition term (<J\) in this case be- 
come larger than the other term. The transition may still 
happen but the complete differential equation has to be 
considered. 

Since an analytic approximation is possible for only 
a limited range of values of (3 in this formulation of 
the problem, the differential equation was also solved 
numerically. We evolve the wavefunction starting from 
p x = — 30v / 2mA to p x = 30\/2mA and the results for 
the tunneling so obtained are shown in Fig. [4] for three 
representative cases of f3 < 0, j3 > and j3 > l/(2m) 
are shown. They show clearly that the oscillations occur 
even outside the limit |2m/?| < 1. Note too that the nu- 
merical data (points in Fig. 4) is in good agreement with 
the WKB result (lines in Fig. 4), except at very small val- 
ues of Eg. The deviation of the numerical solution from 
the WKB estimate at small values of E g (in units of Ao) 
are to be expected since this is precisely the limit where 
the WKB approximation fails. 

Note too that as |/3| is made bigger, the first oscil- 
lation is pushed to ever higher values of E g , and the 
transmission probability at the first antinode becomes 
exponentially smaller. In generic systems, where > 1 
(i.e. the quadratic term in the dispersion dominates over 
the quartic term), the transmission probability at the 
first non-trivial antinode will be extremely small, and 
will likely not be detectable given limited experimental 
resolution. This explains why the common path interfer- 
ence phenomenon has so far been overlooked, despite its 
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universality. However, in bilayer graphene, /3 <C 1 and 
observation of the common path interference may be a 
realistic possibility. 

V. OUTLOOK 

We have shown that common path interference occurs 
generically when the low-order energy dispersion relation 
contains higher order terms in addition to the usual p 2 
term. When the higher order terms come with small 
co-efficients, then the Zcncr oscillations get pushed to 
very large values of the tunneling barrier, such that the 
first antinode involves an exponentially small transmis- 
sion probability. In this situation Zener oscillations may 
not be visible. However, in bilayer (or ABC multilayer 
[H]) graphene, the leading term in the dispersion relation 
is a p 4 term (p 2N ) term. As a result, in the bilayer or 
ABC multilayer, the Zener oscillation effect manifests it- 
self at uniquely small barrier heights, and the first antin- 
ode has unusually high transmission probability. Thus, 
bilayer or ABC multilayer graphene is the ideal place to 
look for oscillations in Zener transmission. 

Other materials with relatively large corrections to 
parabolic spectrum include narrow-gap semiconductors, 
including the recently discovered topological insulators. 
Typical narrow-gap materials are HgCdTe alloys and 
InSb. They are popular infrared photo-detector mate- 
rials. The Kane model is commonly used to describe 
the bands near the T point when the bands are non- 
parabolic. The generic form of bands in this model is 
E = a + bp 2 + \J c + dp 2 , where a, &, c and d are constants 
which can be zero. These bands will always satisfy the 
constraints for common path interference. 

In this paper we considered only the tunneling of elec- 
trons at normal incidence. The full junction IV char- 
acteristic may be obtained by integrating over angles of 
incidence, in the manner of pQ. The resulting IV curve 
are expected to contain regions of negative differential 
resistance, as pointed out in [T], with the negative dif- 
ferential resistance being a consequence of common path 
interference. This may be useful from a device physics 
standpoint. 

Another situation where common path interference ef- 
fects may be important is photon assisted tunneling. In 
the absence of an electric field, the optical absorption 
coefficient goes to zero when photons of the incident 
light have energy below the band-gap. However, in the 
presence of an electric field, absorption of photons can 
take place as shown in Fig. [5| This is known as the 



Franz-Keldysh effect, or photon-assisted tunneling [5]. 
If the bands are simply parabolic, the absorption coef- 
ficient when photon energy hiu < E g is proportional to 

exp(— / ((Eg — hJ)j A ) 3 ^ 2 ), where / is a dimensionless 
constant. When the bands are non-parabolic it is easy 
to see that oscillations should be present in the optical 
absorption below the threshold energy as well. This may 
provide another experimental test of the common path 




FIG. 5: Photon-assisted tunneling of an electron from the 
conduction to the valence band. 



interference phenomenon. 

Thus, we have established that the common path inter- 
ference identified in pQ is in fact a universal phenomenon, 
arising whenever the low energy dispersion relation con- 
tains higher powers of p 2 (but not odd powers of p). 
The phenomenon is easiest to detect when the quadratic 
piece in the dispersion relation is absent (or comes with 
small co-efficient), however, the phenomenon is always 
present. It may be detected in transport experiments, 
or in optical absorption experiments designed to probe 
photon-assisted tunneling. Bilayer (and ABC multilayer) 
graphene are identified as particularly promising materi- 
als for detection of these effects, but other narrow gap 
materials such as several recently discovered topological 
insulators may also offer a playground for exploration of 
these phenomena. 
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